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ABSTRACT 


Buried mines in the beach and surf zone hinder Naval power projection beach 
assault operations and endanger humans. In this research, it is demonstrated that seismic 
interface waves on the surface of a natural beach can be used to detect buried, mine-like 
objects. A seismic sonar system, employing guided interface waves known as Rayleigh 
or Scholte waves has been proposed for detecting shallow, buried ordnance and is the 
subject of an ongoing, ONR-sponsored research program at the Naval Postgraduate 
School. The primary purpose of this research was to develop an improved seismic source 
to evaluate further the concept of using a seismo-acoustic sonar to detect buried ordnance 
in the beach and surf zone. The developed source was based upon a linear magnetic 
force actuator. Additionally, testing procedures were developed for follow-on research. 
Seismic interface waves were generated with two, 25-lb force-controlled linear actuators 
operated as shakers. The waves were measured with a two-element horizontal array of 
three axis seismometers. Buried, mine-like objects, ranging from 71kg to 290kg, and at 
ranges of up to 5 meters were echo-located by visual examination of the raw, recorded 


seismic signal, and by the use of polarization filtering signal processing. 
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I. INTRODUCTION 


Buried mines are a hindrance to Naval beach assault operations. An unknown, 
buried mine, can disrupt, delay, and in some cases, halt Naval operations. In order to 
maintain the relative superiority and capability that U.S. Naval Forces have to project 
power ashore, it is critical that our forces have the ability to counter and neutralize buried 
mines in stride. This capability begins with the detection of buried mines, both in the 
beach and in the surf zone. The proposed concept by the Applied Research Laboratories 
of the University of Texas at Austin (ARL:UT) to use a seismo-acoustic sonar and 
“guided” seismic interface waves to detect shallow, buried ordnance in the beach and surf 
zone offers promise in providing the ability to detect buried mines. This concept will 


further enable our Naval Forces to effectively project power ashore. 


Buried mines are not only a threat to U.S. military forces, they are also a serious 
hazard to civilians and non-combatants. For example, there are more than 100 million 
buried landmines worldwide and more than 150 people are killed or maimed each week. 
In fact, the casualty rate is so high that one in every 470 Angolans and 1 in every 1,000 
Somalis is an amputee. Cambodia alone has an amputee population of more than 30,000. 
(Ref. 1] The concept of a seismo-acoustic sonar that can effectively detect buried mines 
is encouraging and would help reduce the number of landmines and casualties caused by 


buried landmines worldwide. 


The research reported in this thesis is one step towards the development and 
fielding of a seismo-acoustic sonar for detecting buried ordnance. This research is 
presented in seven chapters. The introductory chapter outlines the motivation for 
developing a seismo-acoustic sonar including a brief discussion on continuing plans to 
develop improved Mine Countermeasure (MCM) capabilities for U.S. Naval Forces. 
Also included is a short summary of current and emerging minehunting (MH) methods 
and technologies. The second chapter describes previous research efforts, which helped 
provide guidance for source development and testing procedures. The third chapter 
presents a brief overview on seismic wave theory including a discussion on the 
advantages of using seismic interface waves for detecting buried ordnance. The fourth 


chapter outlines the design and development of the seismo-acoustic sonar and equipment 


used in this experiment. The fifth and sixth chapters describe the field testing procedures 
and experimental results, respectively. The final chapter contains concluding remarks 


and recommendations for further testing and development. 


A. RESEARCH MOTIVATION 


Mine warfare (MIW) has been a part of U.S. warfighting since the American 
Revolution. [Ref. 2] Since that time, U.S. MIW experiences have been infrequent and 
the lessons learned arguably ambiguous. It is partly because of this ambiguity that mine 
countermeasures for U.S. Naval Forces are less than ideal, as evidenced by the damage to 
the USS PRINCETON and USS TRIPOLI by Iraqi mines during DESERT STORM. 
[Ref. 3] 


The reality of an increasing, modern mine threat is not surprising considering that 
most mines are relatively inexpensive and simple to manufacture. The low cost of 
developing an offensive or defensive mining capability combined with its effectiveness 
both tactically and psychologically, make MIW particularly attractive to nations with few 
resources to devote to more expensive, sophisticated weaponry. As a result, mine 
proliferation is becoming a worldwide problem, with more than 49 countries possessing 
mining capabilities, and with 30 of these demonstrating mine production capabilities, and 


20 of these attempting to export their MIW systems. [Ref. 2] 


The increasing mine threat, combined with a downsizing Navy and dwindling 
resources, have driven efforts to develop more flexible and robust minehunting systems 
and capabilities for U.S. Naval Forces. The modemization plan for the U.S. Navy is 
outlined in the United States Naval Mine Warfare Plan [Ref 2] which describes the 
requirement for an organic MCM capability consistent with Navy and Marine Corps 


doctrine, ‘(Maneuver warfare...from the Sea”. 


Power projection operations, including the effective delivery of U.S. Amphibious 
Forces, from the sea, require the avoidance or neutralization of any possible mine threat, 
particularly in the vulnerable regions where amphibious forces transition from shallow 


water, through the surf zone, and onto the beach. Naval mines like the “Manta” bottom 


mine (see Figure 1.1) that damaged the USS PRINCETON typically “scour” below the 
surface and in this process are buried, due to water action. Similarly, beach and surf zone 
mines are purposely hidden or become buried due to erosion and wind action, making 
their detection and neutralization more difficult. These types of buried mines clearly ' 
pose a serious threat to amphibious forces, particularly during an opposed, amphibious 
landing. A flexible, dedicated or organic minehunting capability to detect these buried 
mines in the shallow water surf zone and in the beach 1s needed to effectively neutralize 


the mine threat to U.S. amphibious landing forces. 





Figure 1.1. “Manta” bottom mine similar to the mine that damaged the USS 
PRINCETON during DESERT STORM [Ref. 2] 


Currently, U.S. capabilities to detect and localize buried mines in both the shallow 
water surf zone and in the beach are limited and no one system is optimally capable of 
operating in both environments. The successful development and fielding of a seismo- 


acoustic sonar would provide an additional minehunting capability that could cover the 


critically vulnerable gap from the surf zone to the beach and help neutralize the mine 


threat for our amphibious landing forces. 


B. CURRENT AND EMERGING MINEHUNTING METHODS 


The critically vulnerable regions for amphibious forces as they transition from sea 
to land are the shallow water surf zone and the beach. U.S. amphibious forces are 
vulnerable to buried ordnance in these regions partly because there is no MH system that 
can seamlessly transition from the sea, through the surf zone, and on to the beach while 
continuously hunting for buried mines. While there are a variety of methods and systems 
available to U.S. Naval Forces for hunting mines, no single system is capable of detecting 
buried mines both in the shallow water surf zone and in the beach. Some of the more 


common MH methods are listed below in Tables 1.1 and 1.2 respectively. 









MINEHUNTING METHODS - SHALLOW WATER SURF ZONE 
Navy EOD and SEAL divers 


USMC reconnaissance combat swimmers 





Marine mammals (U.S. Navy dolphins) 
| Towed sweep gear 


Remote Minehunting (RMS) Sonar system 


Airborne Laser Mine Detection System (ALMDS) 


Table 1.1. Current minehunting methods in the shallow water surf zone 





MINEHUNTING METHODS - BEACH 
|Hand-held probes 
‘Manual metal detectors 


‘Ground penetrating radar 


Mechanical breaching equipment 


Explosive breaching nets 


lTrained dogs 





Table 1.2. Current minehunting methods on the beach 


Each of the methods and capabilities listed in Tables 1.1 and 1.2 has distinct 
advantages and disadvantages. For example, manual minehunting methods are generally 
very accurate and have low false alarm rates. However, they are also relatively slow and 
dangerous, as illustrated by the EOD technician probing for buried mines with a hand- 
held probe in Figure 1.2. Additionally, none of the current methods listed in Tables 1.1 


and 1.2 are capable of detecting and locating buried ordnance in the shallow water surf 


zone and in the beach. 
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Figure 1.2. EOD technician with a hand-held probe [Ref. 4] 





There are a number of emerging technologies and systems being considered to 
combat the increasing mine threat to amphibious forces. One of these systems, the 
Coastal Battlefield Reconnaissance & Analysis (COBRA) minefield reconnaissance 
system (see Figure 1.3), an unmanned aerial vehicle (UAV) optical-based sensor system, 
offers promise in that it can detect shallow, underwater mines and is capable of surveying 
large areas. While the COBRA system offers relative improvements in minehunting, it is 


still unable to detect buried ordnance, both underwater and in the beach. 
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Figure 1.3. Overview of the COBRA sonar system [Ref. 5] 


The concept of a seismo-acoustic sonar, as shown in Figure 1.4, offers promise as 
a practical method to detect buried ordnance both in the beach and the seafloor. A 
seismo-acoustic sonar would complement current and emerging minehunting methods 


and technologies, as well as provide an additional sensor for multi-sensor data fusion. 


Like other systems, the seismo-acoustic sonar is highly dependent on the 
propagating characteristics of the medium, which may require the application of 
relatively complex signal processing methods to identify the target signal in the seismic 


sound field. This requirement will be addressed later. 
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Figure 1.4. Design concept for a seismo-acoustic sonar to detect buried mines in the 
shallow water surf zone and in the beach [Ref. 6] 
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Il. PREVIOUS RESEARCH 


There has been a long-standing effort to detect buried mines in the beach and surf 
zone. The Naval Coastal Systems Station in Panama City, Florida has long been working 
on the development of a synthetic aperture sonar that can be “sensor-fused” with 
magnetic sensors for detecting buried mines from the surf zone seaward. Two excellent 
articles on these sensor approaches, one by G.S. Sammelman, et. al. (synthetic aperture 
sonar) and another by T.R. Clem (magnetics), recently appeared in an Office of [Naval 


Research publication Naval Research Review, D. Toderoff, Vol XLIX, b. 3, 1997]. 


Similarly, it is scientifically feasible that seismic interface waves could be used to 
detect and locate shallow, buried, mine-like objects. While there have been many studies 
On seismic waves, there are two recent research studies of particular interest. The first 
study was conducted at the Applied Research Laboratories of the University of Texas at 
Austin (ARL:UT), where field tests were conducted to evaluate the use of seismo- 
acoustic waves to detect a buried mine-like object in a natural beach. The second study, 
conducted at the Naval Postgraduate School, initiated development of a discrete-mode, 
seismo-acoustic source, suitable for detecting buried, mine-like objects in the beach. 


These two research projects laid the foundation for the research reported in this thesis. 


A. ARL:UT (1996) 


Feasibility tests were conducted by ARL:UT to study the possibility of using 
seismic interface waves to detected buried mine-like objects in a natural beach. The 
seismic source consisted of a six-inch by eight-inch exciter foot with protruding nails, 
driven by an electro-mechanical transducer, as shown in Figure 2.1. A three-element 


array of three-axis seismometers was used to receive the seismic signals. [Ref. 7] 


Test results showed that various signal processing techniques were required to 
isolate the target signal in the relatively complex, seismic wavefield. [Ref. 7] These 
signal processing methods included coherent averaging, background subtraction, and 


vector polarization filtering, which will be discussed later. 


The research at ARL:UT validated the proof-of-concept use of seismo-acoustic 
interface waves to detect buried mine-like objects. However, there was observed a 
significant amount of reverberation in the seismic signal. The excessive reverberation 
required the coherent subtraction of the reverberant background to detect the target. 
These results suggested one obvious method of improving the seismo-acoustic sonar: 
reduce the amount of the reverberant background in the seismic, target signals. A logical 
step to reduce the amount of reverberant background noise was to consider a different 
seismo-acoustic source, one that discretely generates only one type of seismic, interface 


Wave. 
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Figure 2.1. ARL:UT seismic target detection beach setup [Ref. 7] 
B. NAVAL POSTGRADUATE SCHOOL (1997) 


Based on the ARL:UT research results, LT Frederick Gaghan studied a discrete- 
mode seismic source designed to selectively excite and generate seismic surface waves. 
[Ref. 8] He developed two discrete-mode seismic sources, the second being a larger 
scale version of the first. The final, field-tested seismic source is shown in Figure 2.2. It 


consists of two, off-the-shelf, 10 pound-force, moving magnet inertial reaction force 


10 


transducers (bass shakers), mounted at 45 degree angles with respect to an aluminum 
base support. Test results showed that, due to medium anisotropy, a drive feedback 
system or different type of seismic source needed to be developed to suppress unwanted 
wave modes. At relatively short ranges, less than about 15 meters (45 feet), it was found 
that multiple body and interface waves were present, resulting 1n a relatively complex 
seismic wave field. At longer ranges, greater than about 25 meters (75 feet), 1t was found 
that the propagating medium naturally filtered the seismic waves so that the wave was 
predominantly a single-mode wave. The most important conclusion relevant to this work 
was that stronger, more linear sources were needed along with refined field testing 


* 


procedures. 





Figure 2.2. NPS discrete-mode seismic sonar design [Ref. 8] 


Based on the conclusions and recommendations from the two previous research 
efforts, the present research focused on the development of a seismic source capable of 


being used to detect and locate buried mine-like objects. Significant efforts were made to 


I] 


refine field-testing procedures and successfully detect a buried, mine-like target to 


support ongoing target strength studies. 
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fil, THEORY OF SEISMIC WAVES 


Seismic waves have been used for geological surveying of petroleum and mineral 
deposits for many years. While there are actually many different types of seismic waves, 
the most common are those generated by an earthquake, which travel through the body of 
the earth and along the earth’s surface. As a result, seismic waves are generally divided 
into two categories, body waves and surface waves. The theory of elasticity predicts two 
basic categories of mechanical body waves for an infinite medium, compressional waves 
and shear waves. In addition to these body waves, there exist many different types of 
surface waves, of which three are of particular interest: the Rayleigh wave, the Scholte 
wave, and the Love wave. Rayleigh waves are surface waves, with both radial (direction 
of propagation) and vertical motion, that exist on the interface boundary between a solid 
and a gas. Scholte waves are similar to Rayleigh waves except that they occur on the 
interface boundary between a solid and a liquid. Love waves are “ducted” waves, 
existing in the layer between three different media or near the interface, where there is an 
increasing shear velocity with depth, and have purely shear motion. This chapter gives a 
brief overview of the basic theory of seismic body waves and surface waves, and 
describes the advantages of using Rayleigh or Scholte waves for detecting buried 


ordnance. 


A. ELASTIC BODY WAVES 


The stress relations for elastic body waves are found by modeling the propagating 


medium as an infinite, homogenous, isotropic and elastic body, and are given by, 


O,, =AG+2ME,, OC, =AG+2pE,, O,, =AG+2QME. 
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where oand € are the stress and strain, respectively, 2 and uz are the Lame constants, and 
@ is the dilatation, or volume strain. The stress and strain indices in (3.1) represent the 


direction of the stress or strain and the plane over which it is applied. PRC I | 
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Neglecting body forces, the equations of motion are given by, 


ae 





d 
p ig =A+u)VO+uV'i, Ge) 
where p is the density of the propagating medium, u represents the particle displacement 


in vector notation, 2 and yw are the Lame constants, and @ is the dilatation or volume 


Strain. [Ref. 9] ° 


The solutions to the elastic wave equations of motion in (3.2) lead to two types of 
mechanical body waves, a compression wave, known as a Primary or P-wave, and a shear 
wave, known as a Secondary or S-wave. The P-wave is the fastest travelling wave, 
producing particle motion that is co-linear with the direction of propagation. P-waves 
can propagate in almost any kind of medium, including rock and water, and travel with a 


longitudinal speed given by 


(3.3) 





where B is the adiabatic bulk modulus, “ is the shear modulus, and pis the density (note 
that the Lame constant w and the shear modulus £ are one and the same) [Ref. 9]. 
Contrary to P-waves, Songucs 46 not propagate in water, due to the fact that liquids 
cannot support static shear stresses. Bulk S-waves produce particle motion perpendicular 


to the direction of propagation and travel with a speed given by 


Sy 


(3.4) 
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where is the shear modulus, and pis the density. [Ref. 9] In general, both P-waves and 
S-waves are produced when an external force disturbs an elastic medium or when an 


clastic surface wave interacts with an interface or a boundary. 
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B. ELASTIC SURFACE WAVES 


While there are actually many different types of surface waves, three specific 
wave types waves were considered in this research: Rayleigh, Scholte, and Love waves. 
As will be shown later, the Love wave was found to result in excessive background 
reverberation, and was not considered further. Rayleigh waves, on the other hand, were 
found to be ideal for detecting shallow, buried ordnance in the beach. While not 
specifically tested, Scholte waves, similar to Rayleigh waves, can also be considered as 


ideal wave types for detecting shallow, buried ordnance in the seafloor and surf zone. 
le Rayleigh Waves 


In 1885, the “great” English scientist, Lord Rayleigh, demonstrated the existence 
of a surface wave that traveled over the interface between a semi-infinite, elastic half 
Space and a gas; hence these waves became known as Rayleigh waves. [Ref. 10] The 
Rayleigh wave equations, obtained by applying the conditions for a stress-free boundary 


to the bulk wave equations in (3.2), are given by 





Q? Q? 
+ =c;V"9, ~ =c;V’y, (3:5) 


where @and ware potential functions representing compression and rotation wave effects 
respectively. [Ref. 9] The finite solutions to (3.5) result in a propagating wave along the 
boundary that is attenuated with increasing depth provided the following condition is met 


Rete |; 


x° —8x* + (24-1677 ? + (l6y? -16)=0. (3.6) 
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The dimensionless variables Kand yin (3.6) are velocity ratios, given by 


ag, 
“" and y= Ll 


c, At+2ue 2-2v 
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where c; is the Rayleigh wavespeed, c; is the shear wavespeed given by (3.4), c; is the 
longitudinal wavespeed, given by (3.3), A and mw are the Lame constants, and v is 
Poisson’s ratio. [Ref. 9] The solution to (3.6) is found for a single, real, positive root on 


the interval x =(0,1), corresponding to a value of Poisson’s ratio on the interval 
v =(0, 0.5). [Ref. 11] For most solid/fluid boundaries, the ratio of Rayleigh wavespeed 


to shear wavespeed is typically about k = 0.9. [Ref. 11] 


As a result of the conditions described by (3.6) and (3.7), Rayleigh waves have 
distinctive features that make them more readily identifiable in a complex, seismo- 
acoustic wavefield. Probably the most obvious feature of Rayleigh waves is the elliptical 
particle motion produced by their passage. This unique characteristic allows Rayleigh 
waves to be distinguished from other wave types by identifying those waves that have a 
90° phase shift (quarter period delay) between their horizontal and vertical components. 
The phase relationship between the horizontal and vertical components resulting in 
elliptical motion is shown in Figure 3.1, which shows an idealized seismic wavetrain 


generated by a single vertical impulse source. [Ref. 12] 


The first two waves (outermost) shown in Figure 3.1a and 3.1b, are the relatively 
fast P-wave and slower S-wave, known also as the “Minor Tremor” in seismic earthquake 
activity. [Ref. 12] The Rayleigh wave closely follows the S-wave and is known as the 
“Major Tremor’. In the Rayleigh wave, the horizontal and vertical components differ in 
phase by about 90°. Near the surface, this results in motion that is counterclockwise 
elliptical, or “retrograde”, with respect to the direction of wave propagation, as shown in 
Figure 3.1c. Note also that, of the three wave types, for a semi-infinite, elastic medium 
with a free boundary, the Rayleigh wave typically has the largest amplitude or 
displacement, particularly the vertical component, which is of importance and will be 


discussed later. 
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Figure 3.1. Seismic wavetrain resulting from a single vertical impulse source [Ref. 12] 


Another important and distinctive feature of Rayleigh waves is that the sense of 
their elliptical motion changes with depth. This is further illustrated by modeling the 
propagating medium, a medium-grain natural beach, as a semi-infinite, uniform (non- 
layered) solid with elastic properties. While this model is not perfectly accurate, it will 
nevertheless illustrate some of the major features of Rayleigh waves. For this present 
project, tests were conducted on a medium-grain, natural beach with a typical Poisson’s 
ratio of about v =~ 0.25, and a measured shear wavespeed of about 90m/s. The seismo- 
acoustic sonar was typically operated at frequencies on the order of about 100Hz. 
Substituting these values into the required condition (3.6), the resulting horizontal and 
vertical Rayleigh wave displacements were determined as a function of depth (see Figure 


3.2) with the MATLAB program in Appendix D23. 


For depths less than about one-tenth of a wavelength (0.2A,), where A, represents 


the Rayleigh wavelength, the horizontal displacement is positive, as shown in Figure 3.2, 
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and the sense of the motion is retrograde with respect to the direction of propagation, as 
seen in Figure 3.lc. At a depth of about 0.2A,, the horizontal motion is zero, resulting in 
purely vertical motion. At depths greater than 0.2A,, the relative horizontal displacement 
changes direction, becoming negative, and the sense of the motion changes to clockwise 
or “prograde” motion with respect to the direction of wave propagation. [Ref. 11] Note 
also that the relative amplitudes of the horizontal and vertical displacements decrease or 
are attenuated with increasing depth into the propagating medium. At a depth of about 
two wavelengths, the amplitude of the displacement is less than 1% of the maximum 
displacement at the surface and for all practical purposes the Rayleigh wave can be 


considered to be localized to within a layer of about two wavelengths deep. 


Plot of Horizontal & Vertical Displacement vs. Depth 
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Figure 3.2. Typical Rayleigh wave displacements in a natural beach as a function of 
depth (see Appendix D23 for the program listing) 
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As shown in Figure 3.1, Rayleigh waves are the largest disturbance in a seismic 
wavetrain generated by a vertical impulse. Additionally, since Rayleigh waves are 
localized near the surface, they attenuate due to geometrical spreading less severely than 
do other wave types. Neglecting absorption and frictional effects, Rayleigh waves 
attenuate approximately as the inverse of the square root of the range (ivr ) from the 
source or scattering center, whereas body waves attenuate approximately as 1/r, due to 
spherical spreading. As a result, Rayleigh waves contain a significant amount of the total 
seismic wave energy, as illustrated in Figure 3.3, which shows the relative energy 
partitions for the Rayleigh, shear, and compression wave components resulting from a 


vertically oscillating circular plate at the surface of a semi-infinite, uniform, elastic 


medium. [Ref. 11] 
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Figure 3.3. Wavefield from a vertically oscillating circular plate at the surface of an 
ideal, semi-infinite, uniform, elastic medium [Ref. 11] 
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Note that Figure 3.3 illustrates many of the Rayleigh wave features that have been 
discussed, namely, that for an idealized, semi-infinite, uniform, elastic medium, Rayleigh 
waves travel along the interface boundary, and contain a significant amount of seismic 
wave energy. Rayleigh waves also result in elliptical particle motion near the surface and 
have a significant vertical displacement component. These distinctive characteristics 
make Rayleigh waves well suited for use in a seismo-acoustic sonar to detect shallow, 


buried ordnance. 
2s Scholte Waves 


Scholte waves are essentially identical to Rayleigh waves except that the interface 
boundary is between a solid and a liquid. Scholte waves contain the same essential 
features as Rayleigh waves, as shown in Figure 3.4. [Ref. 13] Note that Scholte waves in 
the sea floor, similar to Rayleigh waves, are localized near the bottom-water boundary. 
The motion is also depth dependent, changing from retrograde elliptical at the surface to 
purely vertical at a depth of about one-tenth of a wavelength and then to prograde 
elliptical at greater depths. These distinctive characteristics of Scholte waves also make 
them well suited for detecting buried naval ordnance in the seafloor and shallow water 


surf zone sediments. 
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Figure 3.4. Scholte wave characteristics in the seafloor [Ref. 13] 
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IV. DESCRIPTION OF EXPERIMENTAL EQUIPMENT 


The most important part of this research were the research tools, particularly the 
seismo-acoustic sonar. The ‘experimental apparatus is functionally divided into four 
major groups, as illustrated in Figure 4.1: (1) Seismic source(s), (2) Seismic sensor(s), 
(3) Data acquisition, and (4) Signal processing. This chapter describes in further detail 
the system components that were used in this research, including the modifications that 


were made to the seismo-acoustic source, the sensors, and the data acquisition system. 
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Figure 4.1. Simplified representation of the seismo-acoustic sonar system 
A. SEISMIC SOURCE 


The seismo-acoustic sonar source was based on two, commercial, off-the-shelf 
linear actuator systems, made by Aura Systems Inc, (see Appendix Al for specifications). 
[Ref. 14] Depending on the configuration, these actuators are capable of operating as 
either positioners or shakers. For this research, the linear actuators were factory- 
configured to deliver a voltage-controlled force, capable of generating 10lbf-rms and 
25lbf-peak output, with 0.5-inches of total stroke. The actuators have a relatively flat 
(force) frequency response from OHz to 250Hz. The actuators are designed without belts 
or gears so that the only moving part is a moving magnet within a solenoid. The result is 


a simple, sturdy, direct-drive, linear system, without backlash or hysteresis effects. 


I 


The actuators are individually controlled by an Electronic Control Unit (ECU), 
which converts a user-input voltage to a shaft output force. The ECU converts the input 
low-level] analog voltage, to a current that is output to the actuator, typically 1.0-amp 
output per 1.0-volt input (1.0-amp/volt). The actuator output force is linearly 
proportional to the ECU-supplied current. As a result, the force applied at the output shaft 
follows the voltage applied to the ECU, regardless of shaft position. The force constants 
for both actuators used are listed in Appendix A2 and are nominally about 1.29-lbf/amp 


and 1.49-lbf/amp. The linear actuator system configuration is shown in Figure 4.2. 


LINEAR ACTUATOR SYSTEM CONFIGURATION 
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Figure 4.2. The linear actuator system configuration [Ref. 15] 


The actuators can simultaneously generate both DC (static) and AC (dynamic) 
shaft output forces. In addition, there is an internal position sensor, connected to the 
ECU, that allows for automatic centering of the output shaft with up to a 10-pound static 
loading force. The output from the position sensor, which 1s linearly proportional to the 
shaft position, is available via a BNC connector (Position Monitor) on the ECU front 


panel, as shown in Figure 4.3. Tests showed that the position monitor was limited to 
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tracking output shaft movements at frequencies on the order of about 20Hz and below. 
Since the typical operating frequencies used were about 70Hz to 100Hz, a separate 
accelerometer was mounted to the output shaft to track output movements, as shown in 


Figure 4.6. 
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Figure 4.3. Electronic control unit (ECU) front panel 


The ECU requires up to 15 Amps of 11OVAC, 60Hz electrical power for normal 
operation. The actuator was connected to a 19-pin connector (9-pins actually used) on 
the ECU front panel, using a 300-foot underwater cable that was constructed for this 
research (see Appendix B). The underwater cable allowed remote deployment of the 


actuator from the ECU. 
The following CAUTION is of particular importance and should be heeded: 


Do not disconnect the ECU-actuator cable connection without first 
setting the ECU Enable/Disable switch to the “DISABLE” position! 


This lesson was learned first-hand. Disconnecting either the actuator or ECU connection, 
without first turning the enable/disable switch to “DISABLE”, will cause an arc that will 


fuse the cable to the connector. 
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Jl. Actuator Modifications 


Numerous modifications were made to protect the actuators from salt-water 
intrusion and improve their operation as a seismo-acoustic sonar source. Many of the 
modifications were made based on lessons-learned during field tests and by trial and 
error. An external view of the linear actuator, before modification, is shown in Figure 


4.4. 
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Figure 4.4. Schematic of an unmodified linear actuator [Ref. 14] 


One of the first modifications was to attach a nylon cylinder and 5-inch square 
plate assembly to the output shaft, similar to the aluminum rod and plate shown in Figure 
4.7. The actuators were typically placed in an upright position on the beach sand and 
driven with a 5-cycle, 100Hz input signal to the ECU. It was found that for ECU-input 
voltages greater than about 2.5-volts, the actuator-plate assembly literally “hopped” on 
the sand, producing a complex, inconsistent seismic wave at the seismometer, 5-feet 
distant. The cylinder-plate assembly was then buried, in the same sediment conditions, 
about 3-inches below the sand surface. A noticeable improvement was seen in the 
magnitude and the consistency of the seismic wavetrain as seen at the seismometer, 5-feet 
distant. The following week, under similar conditions, a 2-kg inertial mass was fastened 


to the actuator housing (shown in Figure 4.6), which resulted in yet an additional increase 
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in the amplitude of the seismic wave as measured by the seismometer. In this 
configuration, the actuator was operated as an inertial reaction force shaker, with the 
moving mass consisting of the actuator housing plus the added mass fastened to the top 
of the housing. The vertical amplitude of the seismic wave as measured by the 


seismometer, 5-feet distant, is shown in Figure 4.5 for the three modifications. 


Plot of Geophone Output vs. Actuator Input 
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Figure 4.5. Plot of the vertical geophone output voltage vs. ECU input voltage 


As can be seen from Figure 4.5, the addition of an inertial mass combined with 
improved ground coupling by burying the actuator-plate assembly, resulted in an 
approximate gain of 4-dB in maximum achievable seismic wave amplitude. The 
apparent amplitude “rolloff” for the buried plate and mass modifications in Figure 4.5 1s a 
result of the actuator reaching its excursion limit. Note also that the pressure exerted by 


the 25-lbf actuator and 5-inch square plate is on the order of about 1-psi, much less than 
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the liquification pressure for medium-grained sands, which is typically greater than about 


50-psi. [Ref. 16] 


The next series of actuator modifications are shown in Figure 4.6, some of which 
allowed the placement and operation of the actuator inside a watertight, 6-inch diameter 
PVC pipe, shown in Figure 4.7. The modifications were designed to protect the actuator 


from salt-water immersion and allow operation of the actuator system underwater. 


Inertial Mass (2kg) 


ayer 


Re-oriented Connector 90° 


Roller Bearings (4-sets) 


Mechanical Stop 


Rubber Membrane 


Accelerometer 





Figure 4.6. Actuator modifications 


The electrical connector was replaced with an aluminum block that re-oriented the 


connector 90° toward the back of the actuator and allowed placement of the actuator 
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inside the watertight PVC pipe. A watertight, bulkhead connector (see Appendix B3) 
was mated to the PVC housing via a nylon block. Additionally, 4-sets of roller bearings, 
2-rollers each, were attached to the actuator body allowing it to roll relatively freely 
inside the PVC housing. To protect the actuator from internal damage, a hollowed, nylon 
cylinder was mated to the base of the body, over the output shaft, which prevented 


excessive force against internal stops when turning the self-centering actuator off. 


To seal the output shaft end, a nylon ring was mated to the end of the PVC 
housing and 8-holes were tapped into it. Two aluminum plates, with machine screws for 
attaching to the actuator output shaft and ground coupling assembly, were mounted to a 
black, rubber membrane which was then firmly mated via a metal ring to the end of the 
PVC housing. A uniaxial accelerometer (see Appendix A3) was mounted, via a machine 


screw, to the external, aluminum plate, as shown in Figure 4.6. 


To reduce the “head” mass at the actuator output shaft, the nylon cylinder-plate 
assembly, weighting about 1.8-pounds was replaced by a 6-inch long, 0.5-inch diameter 
aluminum shaft and 6-inch aluminum square plate assembly, weighing about 0.8-pounds. 


The fully assembled, watertight, seismo-acoustic sonar source is shown in Figure 4.7. 
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Figure 4.7. The fully assembled watertight seismo-acoustic sonar source 
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om Electronic Control Unit (ECU) Modifications 


After a number of field trips, it was noted that the source continued to “ring” after 
the electrical signal to the ECU was terminated. Laboratory tests showed that the 
actuator continued to oscillate for a few hundred milliseconds at the completion of an 
input pulse train. After consulting with engineers from Aura Systems Inc., it was found 
that the RC time constant and slewing rate of the automatic shaft centering feature had 
been altered by the replacement of the 6-foot, factory-supplied cable with the 300-foot, 
underwater cable. The result was an undesirable oscillation until the actuator settled- 
down to the static, centered position. The solution was to re-adjust the time constant by 
adding a capacitor, in parallel, across the R5 resistor on the main control board. A 
speaker-type push connector was hardwired to the R5 resistor and mounted on the back 
of the ECU as shown in Figure 4.8, enabling external capacitor loading. Trial and error 


lab tests showed that a 2uF capacitor seemed to work best in reducing the oscillations to 


an acceptable level. 
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Figure 4.8. ECU rear cover with an added speaker-type push connector 
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The second modification to the ECU included the addition of an extra pin to the 
19-pin cable connector and tying it to the “spare” BNC connector on the ECU front 
panel, as shown in Figure 4.3. The ECU ts factory-delivered with only 8-pins installed in 
the 19-pin connector, which are all that is necessary to control the actuator. As shown in 
the ECU/Actuator Interconnect Diagram in Appendix A4, the conductors for pins D and 
E on the actuator-side are tied together on the ECU-side before connecting to the ECU. 
After consulting with Aura engineers, pins D and E were tied together at the actuator- 
side, which made one of the 9-conductors in the 300-foot cable available for other uses 


and actually resulted in better grounding and shielding properties. 


Originally, the “extra” conductor, made available from this modification, was 
intended to be used to route the actuator accelerometer signal, which was originally 
installed on the waterproof inside of the PVC housing, back to the ECU and from there to 
the data acquisition computer. Unfortunately, the electrical interference generated by the 
actuator-ECU system added an excessively large 120Hz component to the accelerometer 
signal when using this wire, making it virtually unusable, since typical operating 
frequencies were on the order of 70-100Hz. As a result, a separate coaxial cable was 
used for the accelerometer signal and there is an extra conductor available from the 


actuator to the ECU, which may be used for other purposes, i.e.'a thermocouple. 


It was found that the actuators tended to overheat when constantly operated in the 
vertical position. The overheating was caused by the constant ECU current that was 
required to maintain the statically centered position with the relatively heavy inertial 
mass loaded on top. Therefore the ECU switch was set to “Disable” between 


experimental runs. 
3. Relative Source Gain 


Comparing the vertically operated, linear actuator system with the previous, 
discrete-mode, seismo-acoustic sonar source [Ref. 8], both operating at their maximum 
capability, the actuator produces a much greater seismic interface wave intensity (see 


Figure 4.9). For a source-receiver separation of 70-feet, the actuator source produced a 


Ag) 


vertical seismometer signal of about 60db greater than the discrete-mode source. Based 
on these results, it seems that the linear actuator system looks promising for a seismo- 


acoustic sonar with a reasonable standoff range capability. 
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Figure 4.9. Comparison of the seismic interface waves for the discrete-mode and 
modified linear actuator sonar at a source-receiver distance of 70-feet 


B. SEISMIC SENSORS 


For this research, two three-axis seismometers, consisting of three orthogonally- 
oriented geophone sensors, were used as receiving elements (see Appendix A5 for 
specifications). Each seismic sensor consists of three geophones, mounted in a cylinder 
of urethane (see Figure 4.10) to form a triaxial velocimeter. The geophone converts the 


velocity of the particle motion caused by the propagating wave into an output voltage. 


30 


The output voltage is linearly proportional to the particle velocity, which can then be 
used to further analyze the seismic wavefield, as will be discussed later. The geophones 


have the following the dimensions: height-23 cm, diameter-5 cm. 
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Figure 4.10. Triaxial geophone seismic sensor 
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During field tests there were noticeable frequency components in the seismometer 
signals on the order of about 150Hz to 175Hz, which could not initially be accounted for 
or explained. Similar, unexplainable frequency components were seen during previous 
research studies, using the same geophones. [Ref. 17] Closer examination of the 
geophone specification sheets in Appendix A5 revealed that the geophones have an 
undesirable, “typical spurious frequency” of about 150Hz to 190Hz. Fortunately, for this 
research, typical operating frequencies were below these spurious frequencies, typically 


70Hz to 100Hz. 


The geophone signals are amplified in the deployed unit with a selectable 10dB or 
40dB (typically 40dB selected) amplifier, typically powered by two, 12-volt, 
rechargeable batteries, wired in series to provide +12-volts. During field tests it was 
noticed that the batteries would typically drain to about 10-volts by the end of the day. In 
order to provide more consistent seismic measurements, the batteries were replaced by an 


ac/dc power converter which continuously supplied the required +12-volts. 
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During lab tests it was noticed the actuator-ECU system produced a moderate 
20kHz component in the geophone signal. Using the ac/dc power converter had an added 
benefit in that 1t appeared to reduce the 20kHz interference by about 10%. The most 
effective reduction in the 20kHz noise interference was found by placing the ECUs in a: 
separate equipment rack, effectively isolating the receiving elements from the transmit 


elements. 


€. DATA ACQUISITION 


Data acquisition was accomplished with Signal Processing System’s SPS390 
Signal Analyzer, shown in Figure 4.11. The SPS390 is a 16-bit, 0-40kHz analyzer, 


capable of recording up to eight channels of data. 
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Figure 4.11. Signal Processing Systems (SPS390) signal 
analyzer used for data acquisition 

For this research, the SPS390 was upgraded to a 32-bit operating system and the 
Random Access Memory (RAM) was increased to 1283MB. These upgrades provided 
more a more reliable data acquisition system and allowed data file transfers to another 
computer for signal processing without terminating the SPS390 signal processing 
program. The result was a relatively quick turn-around time (2-3 minutes) from data 
acquisition to signal processing and data analysis. To support follow-on research, the 
configuration, menu Settings, and 8-channel inputs (standardized throughout the research) 


are included in Appendix A6. 
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D. SIGNAL PROCESSING 


The signal processing of the recorded seismic data during this research has been 
an evolutionary process. At the start of this research, seismic data were recorded and 
stored in the field, then processed in the lab 2 to 3 days later. This method was found to 
be excessively inefficient, considering that it would take a week or more to analyze the 
data. As a result, equipment and signal processing capabilities were developed to allow 
near real-time (2-3 minutes delayed) signal analysis in the field. This capability allowed 
an on-site assessment of the effectiveness of sonar placement and test configurations, 


significantly improving testing procedures and results. 


Twenty two signal processing subroutines (see Appendix D) were written using 
Matlab’s Signal Processing Toolbox. [Ref. 18] The data analysis and signal processing 


methods that were developed are summarized in 4 basic steps, shown in Table 4.1. 


LL = 


OPERATION 


SPS390 Record accelerometer and seismic geophone signals 


SPS390-to-Portable | Download the recorded data (via a null-model cable) 


Portable ; Signal Analysis (MATLAB Subroutines) 





‘i Portable Convert the data to Matlab 4.1 format | 


Table 4.1. Summary of data acquisition and signal processing steps 


The MATLAB subroutines, listed in Appendix D, are sufficiently commented and 
do not require further explanation at this point. One particular signal processing method, 
vector polarization filtering, deserves further discussion. As discussed in Chapter II, 
Rayleigh waves, by their very nature, have radial and vertical components that differ in 
phase, resulting in polarized elliptical particle motion. This feature provides a method 
for discriminating between different wave types to identify those interface waves with a 
relative phase difference between the radial and vertical components. This method is 


known as “Vector Polarization Filtering”. [Ref. 7] 
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In this research, polarization filtering is the process of extracting from a radial- 
vertical signal pair, the component of the crossed signal power that comes from the 
elliptical particle motion of the interface wave. [Ref. 7] This is accomplished by 
applying the Hilbert Transform to the radial-vertical signal pair (see Ref. 7 for further 


discussion) and calculating the time domain, crossed power function as given by, 
Fy = Nimes *Vaitbert (4.1) 


where P, is the crossed power function, r,.,.,is the complex conjugate of the hilbert 
transform of the radial component of motion, and y,,,.,, 1S the hilbert transform of the 


vertical component of motion, consistent with the notation used in Ref. 7. The in-phase 


radial and vertical components are contained in the real parts of P, while the out-of- 


phase components are contained in the imaginary parts. [Ref. 7] As a result, the 


imaginary part of P, is approximately proportional to the power of the Rayleigh wave 


and deflects positively or negatively for polarized signals, depending on whether the 
motion is prograde or retrograde, respectively. [Ref. 7] This technique is used to reject 


other signals, particularly the bulk P and S waves, and will be illustrated later. 
Fr. TRANSPORTATION 


During this résearch a John Deere, 4x6 “TrailGator” (see Appendix A7) was 
acquired to move and position naval mine-like targets in the beach, see Figure 4.12. 
Additionally, since part of the focus of this research was to refine testing procedures, 
considerable efforts were made to minimize and simplify equipment requirements. As a 
result, the transmit, receive, and signal processing equipment components described 
previously were condensed down to two equipment racks that were transported and 
staged from the bed of the TrailGator, as shown in Figure 4.13. The TrailGator proved to 
be an invaluable addition to the research project and presents opportunities for expanded 
development (i.e. deploying a combination seismo-acoustic sonar and sensor array), 


which will be discussed later. 
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Figure 4.12. Front view of the TrailGator and beach equipment 
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Figure 4.13. Transmit and receive equipment racks in the TrailGator bed 
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V. TESTING PROCEDURES 


One of the major goals of this research was to refine field-testing procedures for 
the seismo-acoustic sonar development. As a result, a significant amount of time and 
effort, more than seven months, was put into developing streamlined and efficient beach 
setup and testing procedures. Rather than describing the evolution of the developing test 
procedures, this chapter will outline the procedures that were developed and successfully 
used at the later latter part of this research. Before describing the testing procedures, this 
chapter begins with a brief description of the beach test site, characterization of the beach 


sand, and a spectral analysis of the ambient ground vibration. 


A. BEACH TEST SITE 


The beach test site was located on the Monterey Bay coastline, adjacent to the 


Naval Postgraduate School, at the NPS Beach Research Facilities. See Figure 5.1. 
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Figure 5.1. Beach and underwater test sites 
Map copyrighted 1994 by the California State Automobile Association. Reproduced by permission. 


The beach site proved to be well suited for this research, with easy access to the 
beach via a 20-foot wide, dirt road. The width of the beach and available test area varied 


with the tides, ranging from a width of about 120-feet at low tide to less than 10-feet at 


highsiide. See Figure 5.2. 
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Figure 5.2. Beach test site at high tide 


By. SEDIMENT CHARACTERIZATION 


In general, beach sand 1s composed of a variety of elements, including minerals, 
rocks, organic and inorganic compounds, water, and air. The relative amounts of these 
different elements vary, depending on the geographical location and the depth of the 
sediment. [Ref. 8] In addition to its composition, sand is characterized by a variety of 
features, including grain size, density, and moisture content. Each of these characteristics 
affects the propagation of elastic wave energy along the beach-air interface. While the 


relationship between sand characteristics and interface wave propagation are relatively 


38 


complex, it is generally accepted that interface wave energy is less attenuated in high 


moisture content and densely compacted soil environments. [Ref. 16] 


Although the Rayleigh wave equations, for a uniform, semi-infinite, elastic half 
space (3.5), are frequency-independent (dispersionless), real beach sediments are not 
uniform, meaning that they are not necessarily dispersionless. Hence the speed of 
Rayleigh waves may depend on the frequency. In a natural beach or seafloor sediment, 
the overburden weight of the surface sediments causes an increase in shear modulus with 
increasing depth into the medium. Consequently, longer wavelength, lower frequency 
acoustic waves propagate deeper into the sediment, where the shear velocity is greater. 
AS a result, long wavelength, low frequency signals travel faster than short wavelength, 
high frequency signals. At relatively higher frequencies, greater than 100Hz, there is little 
dispersion because the wave remains near the surface, where the sediment is relatively 
uniform. To complicate matters further, a core sample of the beach sediment revealed, 
unsurprisingly, that the moisture content of the sand varied with depth. At a depth of 
about 1-meter, it was found that the sand was completely saturated, corresponding to the 


depth of the tide water table at that time. 


C. AMBIENT NOISE SPECTRA 


Using the geophones as receivers, a measurement of the ambient ground vibration 
at the beach test site was made. The measurement allowed an assessment of the influence 
that the ambient noise would have on identifying the target signal, and provided 
qualitative information on which operational frequencies to avoid. See Figure 5.3. The 
test setup is outlined in Appendix Cl] and the MATLAB program used for the spectral 
analysis is outlined in Appendix D16. 


Note that the spectral content of the ambient background noise is greatest from 
about 10Hz to 40Hz, consistent with previous measurements made at this site. [Ref. 3] 
The predominant source of the ambient noise is undoubtedly the breaking of waves 
similar to those shown in Figure 5.2. The measurement was useful in providing guidance 


on selecting the high-pass filtering frequency and allowed a relative sense of confidence 
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that the sonar operational frequencies, typically 70Hz to 100Hz, are not within the strong 


ambient noise band. 
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Figure 5.3. Omnidirectional velocity power spectral density for the ambient ground 
vibrations at the beach test site (see Appendix Cl for the test setup) 


SOURCE EMPLOYMENT 


Two seismic sources were developed to independently generate horizontal and 


vertical motion. The versatility of the sources allowed for creative testing of many 


different actuator orientations and operating combinations. More than 100 different tests 


were conducted over 21 days of field testing, each with different drive signals, actuator 


orientations, geophone locations and target positions. The seismic interface waves 


resulting from these many tests were closely studied and evaluated for further use in 


buried target detection. 
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The placement and operation of the actuators can be categorized into two different 
configurations. In the first configuration, the actuators are operated either horizontally or 
vertically, while in the second, they are operated simultaneously, one horizontal and one 
vertical, with a 90° relative phase difference between their drive signals. With the reader — 
as the receiver, Figure 5.4 shows a typical vertical and horizontal actuator placement, 


resulting in vertical and shear wave excitation, respectively. 


—- 


Figure 5.4. Typical actuator placement for vertical and shear wave excitation. 
1. Placement and Orientation 


Different actuator orientations and operating combinations resulted in the 
generation of various seismic interface waves, each with distinctive features and 
characteristics. These waves were evaluated for their suitability for use in detecting 
buried targets. Based on previous research [Refs. 7 & 8], the evaluation centered on 
finding seismic waves that had both a single-mode polarization and minimal 


reverberation. 


The first, and simplest, actuator placement that was tested was to orient the source 
vertically, as shown in Figure 5.4. The test setup is shown in Appendix C2, which shows 
two actuators separated by 46-cm (18-inches) and 1.8-meters (6-feet) from the receiver. 
The actuators are driven with a single cycle, 10-volt, 90Hz signal. The geophone signals 
for the relative x, y, and z-motion are shown in Figure 5.5, which was generated using the 
MATLAB program in Appendix D4. The vertical line at about 0.01-seconds indicates 


the start time for the input actuator pulse. Note that the x-radial and z-vertical motions 
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are the dominant wave components and that the reverberant wavefield has decayed to less 


than 5% of the maximum amplitude at about 0.1-seconds. 
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Figure 5.5. Relative seismic motion (velocity) from a vertically-oriented source 


As discussed in Chapter II, computation of the crossed power function (4.1) is a 
practical means of separating the Rayleigh wave signal and rejecting the bulk P and S 
waves. A plot of the imaginary component of the crossed power function provides an 
indication of the polarization, or direction, of the elliptical wave motion, prograde or 
retrograde. With the seismometer reference post (see Figure 4.10) pointing in the 
direction of the incident wave, a positive deflection of the crossed power from (4.1) 
indicates that the motion of the Rayleigh wave is prograde, whereas a negative deflection 
indicates a retrograde wave. The observed imaginary crossed power for a vertically- 


oriented source is shown in Figure 5.6, which'was generated using the MATLAB 
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program in Appendix D12. Note that the deflection is negative, corresponding to an 
interface wave with prograde motion for a geophone that is pointed away from the 


incident wave source as shown in the test setup in Appendix C2. 


Similar tests were conducted for a horizontally-oriented actuator first pointing 
directly at the receiver at a 0° angle, resulting in a P-wave, and then transversely at a 90° 
angle, resulting in a Love wave. In both cases it was found that reverberations were 
excessively large and that there were at least two excited modes, as the polarization of the 
elliptical particle motion was observed to change from one sense to the other within the 
pulse train. These two actuator orientations were not considered further for use in 


detecting buried targets. 
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Ficure 5.6. A plot of the imaginary crossed power for a vertically-oriented source 
with the receiver directed away from the source 
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2: Phasing 


In the second type of configuration, the actuators were oriented side-by-side, one 
horizontal and the other vertical. The actuators were driven with a relative 90° phase: 
difference to simulate and excite either retrograde or prograde Rayleigh wave motion. 
The results were ambiguous; in all of the tests the polarization of the elliptical particle 
motion changed during the pulse wavetrain. Similar phenomena had also been observed 
during previous studies [Ref. 8]. It is unclear whether there were actually two different 
modes, One prograde and one retrograde wave, or whether the same wave was being 
reflected from the saturated layer about 1-meter in depth. It seems unlikely that the 
change in polarization was due to reflection from the underlying layer, as the wave 
polarization for the vertically oriented source should have been just as likely to generate 


waves that were reflected from the underlying layer. 


The test configuration for generating retrograde and prograde motion with the 
actuators is shown in Appendix C3. The test configuration shows two actuators, one 
horizontal and the other vertical, separated by 46-cm (18-inches) and driven with a two 
cycle, 6-volt, 100Hz signal with a relative phase difference between the two, for 
simulating retrograde or prograde elliptical motion. Figure 5.7 shows the crossed power 
calculated from (4.1), which indicates the polarization or direction of motion for a 
retrograde input and a Source-receiver distance of 5-feet (1.5-meters) and 10-feet (3.0- 
meters). Note that the polarization initially starts out as prograde at 5-feet (1.5-meters) 


and changes to a mixed retrograde/prograde wave at 10-feet (3.0-meters). 


Using the same geometry and configuration as for the retrograde test, the phase of 
the input signal was changed to simulate and excite prograde motion. The resulting 
motion is indicated in the plot of the imaginary crossed power in Figure 5.8. Note that 
the polarization is a mixed retrograde-prograde motion at 5-feet and changes to purely 
prograde at 10-feet. Although there appears to be a single-mode, prograde polarization at 
10-feet, similar to the polarization for the vertically oriented source, this type of wave 
was not used because of the desire for a short-range sonar system and the fact that the 


mechanism for mixed polarization was not fully understood. 


44 


Source-Receiver Distance (5 feet) 


anh 
—~ oO 


Amplitude (volts) 
o) 





0 0.05 0.1 O15 0.2 0:25 
Time (sec) 


‘Source-Receiver Distance (10 feet) 


= 
on 


Amplitude (volts) 
© 





0 0.05 0.1 0.15 0.2 0.25 
Time (sec) 


Figure 5.7. The imaginary crossed power resulting from a retrograde-phased 
actuator pair with the receiver oriented towards the source 


In all tests of this type, orienting the two actuators, one horizontal and one 
vertical, and driving them with a relative phase lag resulted in a combination of modes or 
wave types that included both retrograde and prograde motion. After extensive testing, it 
was determined that the vertically oriented source resulted in a simpler, single-mode 
wave that was more suitable for detecting buried targets. Compared with other wave 
types, the vertical source produced waves that had relatively less reverberation and a 
consistent, single mode or polarization. The following testing procedures are based on 


operating the seismo-acoustic sonar as a vertical source. 
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Figure 5.8. The imaginary crossed power resulting from a prograde-phased 
actuator pair with the receiver oriented towards the source 


E. TESTING CHRONOLOGY 


The bulk of the equipment was transported from NPS to the beach site in a large, 
unfurnished mobile home. The mobile home was stationed on the bluff overlooking the 
beach and served as a Staging area as well as supplying electrical power for the beach 
equipment. The equipment, including mine-like targets, was transported from the bluff to 
the beach using the TrailGator 4x6 utility vehicle. Once on the beach, the equipment was 
unloaded and set up according to pre-planned test configurations as shown in Figure 5.9. 
After the initial equipment setup, optimization measurements were made to find the best 


operating frequency and the associated wavespeed at that frequency. 
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Figure 5.9. Equipment setup at the NPS beach research site 
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1. Frequency 


Based on previous research, the driving frequency for the source was initially set 
at 100Hz. Typical velocities for this frequency were found to be on the order of about 80 
m/s +10, resulting in a wavelength of about 0.8 m. After many unsuccessful target 
detection attempts at this frequency, consideration was given to testing other frequencies. 
In general, it was found that frequencies below 100Hz worked best, consistent with 
previous “ground truthing’” measurements made to determine the frequencies and 


wavespeeds of the dominant interface wave modes. [Ref. 8] 


Generally, it was found that at frequencies of about 55Hz and below, the motion 
of the drive plate, as measured by the accelerometer, became noticeably different from 
the drive waveform. The cause of this apparent non-linearity is not known. At higher 
frequencies, greater than about 150Hz, the interface wave was found to quickly attenuate 
with range. Based on trial and error, various frequencies, starting from about 60Hz and 
increasing in 10Hz increments, were tested on-site to determine the frequency that 


resulted in the wave with the largest vertical component. 


The evolving waveform as a function of distance for two, different operating 
frequencies, 100Hz and 70Hz, are shown in Figures 5.10 and 5.11, respectively. Note 
that the wavespeed for the 70Hz wave is slightly higher, as it penetrates.deeper into the 
sediment, where the shear velocity is higher. More importantly, the 70Hz wave is 
attenuated less than the 100Hz wave, making 70Hz the optimal source drive frequency 
for the sediment on this particular day. In general, the optimal frequency was determined 
by placing the geophones at a fixed distance from the source, typically 40 feet (see 
Appendix C4), and adjusting the frequency to find the wave with the largest vertical 


amplitude. 
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Figure 5.10. Surface wave evolution for a 100Hz input drive signal 
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Figure 5.11. Surface wave evolution for a 70Hz input drive signal 
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2. Wavespeed 


After determining the optimal operating frequency (test setup in Appendix C4), the 
wavespeed was measured by cross correlating the vertical geophone signal with the 
source accelerometer signal (see the MATLAB program listing in Appendix D7). The 
cross correlation provided the time delay for maximum correlation between the two 
signals, providing the time that it took for the wave to travel from the source to the 
receiver. Based on the known source-receiver distance and the measured time delay, the 
wavespeed was quickly determined, within 2-3 minutes. The measured wavespeed was 
then used for echo-ranging in subsequent background measurement and target detection 


tests. 


3. Background Data 


After determining the optimal drive frequency and its associated wavespeed, the 
geophones were repositioned at an appropriate distance from the source, typically 6 to 10 
feet (1.8 to 3.0 meters), as shown in Appendix C5. A series of 10 background (i.e. no 
target) measurements were made using a single cycle drive cycle at the optimal 
frequency, and the geophone signals coherently averaged. This averaged “background” 
data was then analyzed using the vector polarization filtering algorithm in Appendix D12 
to determine the time that the wavefield, as seen at the geophones, had decayed back to a 
nearly quiescent state. Based on the known wavespeed, the minimum range at which a 
buried target might be detected was estimated. Using this information, the target was 


buried at an appropriate distance from the receiving geophones. 
4. Target Data 


Two differently shaped and sized targets were tested in this research. The mass of 
these two targets was varied, and will be discussed later. The target data collection test 


phase was conducted by burying a target, typically 5 cm below the surface, at an 
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appropriate distance, and collecting 10 data sets, using the same drive conditions as for 
the background data. As necessary, the data were coherently averaged and vector 
polarization filtered (see Appendices D2 and D12 for MATLAB subroutines) for study 


and target detection analysis. 
5. Signal Processing 


Numerous signal processing routines, more than 23 total, were developed, using 
MATLAB’s signal processing toolbox. [Ref. 18] These subroutines were invaluable in 
allowing near real-time analysis of the measured seismic signals. As previously 
discussed, the background and target signals were coherently averaged and polarization 
filtered. Also available for target signal analysis was a coherent subtraction subroutine 
(see Appendix D19), which coherently subtracted the background, reverberation field 
data from the target signal. This capability allowed a comparison of before and after 
seismic wavefield characteristics to determine. and assess the detection of the target 
signal. Also available was an incoherent subtraction subroutine (see Appendix D18) 
which subtracted the background imaginary power from the target imaginary power. 
Similar to the coherent background subtraction, this method offered an additional 


capability to assess detection of the target signal. 


oy | 








VI. TEST RESULTS (VERTICAL SOURCE) 


Summarized in this chapter are some of the important test results from a vertically 
oriented seismic source. As previously mentioned, because of the versatility of the | 
seismo-acoustic actuator source, it was possible to conduct numerous tests studying 
seismic interface wave generation and propagation. The results presented in this chapter 
include wave attenuation, source velocity, energy, and displacement measurements, 


target detection, and underwater testing. 
a: ATTENUATION 


Although it was not studied in depth, limited attenuation measurements were 
made to help evaluate the types of seismic waves that were being generated. Based on 
the theoretical results presented in Chapter III, neglecting absorption, Rayleigh interface 
waves can be expected to attenuate due to cylindrical spreading. As a result, the 
acoustical power of the cylindrically spreading interface wave can be expected to 
decrease as the inverse (1/r) of the radial distance from the source. Since the complex 
power (4.1) gives the relative power of the seismic pulse [Ref. 7], determining the 
attenuation of the seismic wave is relatively straightforward. Using the 70Hz seismic 
source and evolving seismic pulse shown in Figure 5.11, the cross power function was 
calculated and the maximum observed values were plotted as a function of range in 
Figure 6.1. Curves for cylindrical (1/r) and spherical spreading (1/r*) are also included 


for comparison. 


The measured attenuation curve in Figure 6.1 indicates that the seismic wave is 
attenuated less severely than for spherical spreading. In addition to spreading effects, 
absorption from the friction of the sand grains rubbing together results in wave 
attenuation. The results shown in Figure 6.1 qualitatively suggest that the vertical source, 
operating at 70Hz, is generating an interface wave that attenuates due to cylindrical 


spreading. 


53 


10° 


Measured 


R* Dependence 
R" Dependence 


10 


fe 


= 
oO 
© 


Cylindrical 
Spreading 


RMS Power (P 


10 


Spherical a 
Spreading 





10, 
4 6 8 10 12 14 16 18 20 22 24 
Range (m) 


Figure 6.1. Seismic wave attenuation using two vertical actuators operating at 70Hz. 
B. VELOCITY AND DISPLACEMENT 


The velocities and displacements of the motion produced by the vertically 
oriented seismic source were determined for the test setup shown in Appendix C2. The 
recorded geophone signals are shown in Figure 5.5. The geophone velocities are given 
directly by the output signal and the source velocity was found by integrating the actuator 
accelerometer signal. The displacements were obtained by integrating the velocities. 
The relative velocities and displacements for the seismic source and geophone (computed 
sing the MATLAB subroutines in Appendices D20 and D8, respectively), 6-feet (1.8- 


meters) distant, are shown in Figures 6.2 and 6.3, respectively. 
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Figure 6.2. Relative velocities of the linear actuator shaft and geophones 
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Figure 6.3. Relative displacements of the linear actuator shaft and geophones 
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The average power transmitted from the actuator to the sediment was estimated 
by calculating the product of the actuator force times the drive shaft velocity. It was 
determined that approximately 0.6-watts of seismic power was generated at the drive 


plate of a single, linear actuator system. 
C. TARGET DETECTION 


After months of testing and equipment modifications, two buried mine-like 
targets were unambiguously detected using the linear actuator seismo-acoustic sonar 
system. The first target was a 156-lb (71-kg) helium gas tank, shown in Figure 6.4. The 
gas tank was modified so that the top could be removed to allow for the placement of 26- 


lb (11.8-kg) lead blocks inside, thus increasing its mass. 





Figure 6.4. Buried helium gas tank target 


The second target, shown in Figure 6.5, was a |.5-feet (0.68-meter) diameter, 1|- 
foot (0.45-meter) height gunpowder keg weighing about 16-Ilbs (7-kg) empty. The 


removable lid allowed placement of 26-lb blocks inside, shown in Figure 6.5. 
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Figure 6.5. Buried gunpowder keg target with 25-lb lead blocks alongside 


A test was conducted to detect the buried helium tank in Figure 6.4. The test 
setup, outlined in Appendix C5, consisted of two linear actuators, oriented vertically, and 
driven with a single cycle pulse at 70Hz. The geophones are 5-feet (1.5-meters) from the 
source and the target is buried 11-feet (3.35-meters) from the source. The wavespeed 
was determined to be about 290-ft/s (88-m/s). The background and target signals, using 
10 coherently averaged records (subroutine in Appendix D2), for the source, and the 
relative the radial, shear, and vertical components of the geophone motion are shown in 
Figure 6.6. Note that the target signal, which should be seen at about 0.08 seconds, is not 


readily identifiable in the raw, recorded data. 


The crossed power for the background data, target data, and target data with 
background coherently subtracted out was calculated using (4.1) and the results are 
presented in Figure 6.7. The measured wavespeed of the day (67m/s) was used to 


transform the time-based record to a range-based record. 
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Figure 6.6. Relative source and geophone motions for the background and target data. 
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Figure 6.7. Imaginary crossed power for the background and target data. 
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Relative Motion (volts) 


Note that the imaginary crossed power gives a clear indication of the presence of 
the buried tank at about 17-feet (5.2-meters), as expected for a round trip range of 11-feet 
(3.3-meters) to the target and 6-feet (1.8-meters) back to the geophone. The negative 
deflection of the imaginary power indicates that the motion of the target echo is 


retrograde. 


Similar tests were conducted two weeks later. The tests used both the helium tank 
and the gunpowder keg, which were buried a few days prior to the tests. In these tests, 
the geophones were positioned at 10-feet (3.0-meters) and the targets at 16-feet (4.9- 
meters). The test setup is outlined in Appendix C6. The mass of the targets was also 
increased in 104-Ilb (47-kg) increments using the lead blocks shown in Figure 6.5. The 
RMS values of the relative motion for the buried, 572lb gas tank are shown in Figure 6.8. 


Note the apparent evidence of the target at about 23-feet (7-meters). 
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Figure 6.8. Relative RMS source and geophone motions for the buried 572]b gas tank. 
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The crossed power for the gas tank with increasing mass loading was calculated 
using (4.1) and the results are presented in Figure 6.9. Note that the imaginary crossed 
power gives a clear indication of the presence of the buried tank at about 23-feet (7- 
meters) as expected for a round trip range of 16-feet (3.3-meters) to the target and 6-feet 
(1.8-meters) back to the geophone. The relative target strength appears to be related to 
the target mass, as seen by the increasing imaginary power amplitude with increasing 
target mass. The positive deflection of the crossed power indicates that the motion of the 
target signal is prograde, which is surprising. It should be noted that, both targets had 
scoured themselves to a depth of about 0.5-feet (0.15-meters) due to water action from 
relatively high tides. It may be that this results in target scattered signal motion that is 


prograde. 
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Figure 6.9. Imaginary crossed power for the helium gas tank target 


Similar results were seen for the gunpowder keg, which was left buried in the 


beach for 4-days after the helium gas tank measurements. The crossed power for the 
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gunpowder keg with increasing mass loading was calculated and the results are presented 


in Figure 6.10. The results are similar to the buried gas tank, which was detected four 


days earlier. The presence of the keg is clearly indicated by the positive deflection of the 


crossed power at a range of about 24 feet (7.3-meters). As with the gas tank, the relative 


target strength increases with increasing target mass. 
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Figure 6.10. Imaginary crossed power for the gunpowder keg target 


Relatively unambiguous target detection using the vertically oriented seismo- 


acoustic source was seen for both the helium gas tank and the gunpowder keg. 


D. 


UNDERWATER TEST 


Near the completion of this study, a single field trip was made to test the seismo- 


acoustic sonar system underwater. Unfortunately, high tides and storm surge prevented 


6] 


successful deployment and underwater testing. The lessons learned from the test are 


included to help support further testing in this area. 


The underwater testing initially started at the beach site. It was found that the 
300-feet underwater actuator and geophone cables were not long enough to allow 
deployment just beyond the surf zone. Although not available for that particular day, a 


small boat would allow deployment of the system near the surf zone, from the water side. 


Since deployment from the beach was not possible, the equipment was loaded 
into the mobile home and transported to the municipal wharf at the Monterey Marina. 
See Figure 5.1. Unfortunately, the relatively strong storm surge did not allow underwater 
placement of the actuators or geophones. Typically, this site would be ideal for 
underwater tests, as it is known for its relatively good visibility and flat, sandy beach 
bottom. The location is also ideal because it allows use of the mobile home to stage the 


electrical equipment and has easy beach access for water entry, as shown in Figure 6.11. 
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Figure 6.11. Underwater test site showing easy water access (the author and his spouse) 


VII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


A uniaxial seismo-acoustic sonar source based on a linear force actuator was 
developed for this research. This source resulted in an impressive 60db improvement 
over the previous seismic source that was developed at NPS. [Ref. 8] Previous efforts at 
ARL:UT resulted in the detection of a buried object at a range of 3-meters after 
coherently subtracting the reverberant background data from the target data. [Ref. 7] 
Clearly it is not practical to rely on background measurements to detect buried mines in 
an operational environment. This source detected buried mine-like targets at ranges up to 


5-meters and did not require coherent background subtraction. 


The linear force actuators worked exceptionally well to generate seismic interface 
waves with minimal reverberation. The complexity of the seismic wavefield was reduced 
and buried target detection was enhanced. The results of this thesis support the use of 


seismic interface waves to detect ordnance. 


B. RECOMMENDATIONS 


While the linear force actuator source resulted in successful target detection, it is 
limited in the amount of force that it can apply to the sediment (25-lb peak force). 
Therefore, the intensity of the seismic interface wave and the effective operating range of 
seismo-acoustic sonar system is limited by the force capability of the actuator-based 
system. To increase the intensity and range of the system, a single, waterproof actuator, 


capable of 100+lb peak force should be acquired and used in subsequent studies. 


In addition to a more intense source, an array of receivers should be developed to 
determine the azimuth or bearing to the target. The combination of an intense seismic 
source combined with an array of geophone sensors will allow for bearing and range to 


target measurements at greater standoff ranges. 
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Along with the development of a seismic receiver array, a MATLAB or other 
appropriate subroutine should be written to make practical use of the range and azimuth 
data provided by the array. Conceivably, within 2-3 minutes after recording the data, a 
polar plot display, similar to the display presented in Figure 7.1, would allow for a near 
real-time assessment of buried mine-like objects in the seismic wavefield. The display 


could also be updated as more data was collected and processed. 


Figure 7.1. Polar plot display showing bearing and range to a buried target [Ref. 7] 
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The SPS390 Data Acquisition and Signal Processing System has the capability to 

store and average data based on an external trigger signal. To support increased and 
more effective data collection, the averaging feature available with the SPS390 should be 
used. The SPS390 owners manual [Ref. 19] contains an “Averager” conversion program 
that would allow for the conversion of the .avg file format into an ASCII file format 


compatible for analysis using MATLAB’s Signal Processing tools. 
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Finally, underwater tests should be initiated as soon as possible. In theory, 
underwater Scholte waves are nearly identical to Rayleigh waves. The seismo-acoustic 
source developed in this research should work equally as well underwater as it did on the 
beach. However, this needs confirmation. In fact, since naval mines are typically more 
massive (500lbs to 2000lbs) than land-based mines, it should relatively easier to detect 
them using a seismic-based sonar. Conceivably, an autonomous underwater vehicle 
could be used as the supporting platform for a seismo-acoustic sonar system to detect 


buried naval mines in the shallow water near the surf zone and seaward. 
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APPENDIX A. EQUIPMENT SPECIFICATIONS AND SPS390 SETTINGS 


This appendix contains the equipment specification sheets and SPS390 data 


acquisition settings that were used during this research. 


Appendix Al — Aura linear actuator and ECU specifications 
Appendix A2 — Linear actuator force constants 

Appendix A3 — Actuator accelerometer specifications 
Appendix A4 — ECU/Actuator interconnect diagram 
Appendix A5 — Geophone specifications 

Appendix A6 — SPS390 setup and menu settings 

Appendix A7 — John Deere TrailGator specifications 


67 


APPENDIX Al. LINEAR ACTUATOR AND ECU SPECIFICATIONS 
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APPENDIX A2. LINEAR ACTUATOR FORCE CONSTANTS 


FORCE IN POUNDS 





CURRENT IN AMPS 


FORCE CONSTANT: 4 to 16 + @ = 1.49 
COMMENT: UUT MOUNTED VERT 1.7H ADDED TO Y AXIS, MID-STROKE 


FILE: AL@1F161.DAT SETUP: A12S-FC.DSU 10-24-1997 11:83 


A101A ACTUATOR FORCE CONSTANT SN 162 


FORCE IN POUNDS 





CURRENT IN AMPS 


FORCE CONSTANT: 4 to 16 + 0 = 1.2 
COMMENT: UUT MOUNTED VERT 1.7 ADDED TO Y AXIS, MID-STROKE 
FILE: AN@IFIGe.DAT SBWUR2"AleS-ec DSU 10-24-1397 10:54 
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APPENDIX A3. ACTUATOR ACCELEROMETER SPECIFICATIONS 


with built-in electronics.§ | 8.5 Eh 


Dimensions shown in incnes except where noted. 


10-32 Thd Coaxial 
5-44 Thd Coax. Connector 
5/16 Connector j = 
Hex aap be 
I t .65 
33 | 
al a i l 

.10 Dia .10 Dia Typ 

5-40 Thd 5-40 Thd Hex 10-32 Th 

Integral Mig Stud Integral Mig Stud Mtg Hole 


NO 









351B11 351B14 
(351812) Top connector mode! numbers listed in parentheses 















10-32 Thd Coax. 
Connector 
.50 


351B03 


(351804) 





HIGH FREQUENCY / LOW MASS 


te ater sa5% 
; Bre ee By eee EER at eaens (se 


| English Selection Guide 


~ .75 Dia a 






me 75 Dia 7 







oe : 10-32 Thd fhe T— 10-32 Thd 
lea oi | Coaxial Coaxial 
| Connectorcto Connector 
85 85 





7” = | 
34 =_ 25 aa 
Hex 40.32 Thd Hex 140-32 Tha 
Mtg Hole Mtg Hole 
351831 351841 
351832) (351B42) 





GENERAL 


HIGH SENSITIVITY 





MODEL NUMBER [unit [asian [| —astere | —ssveca «|= a51831 «| ~—sastBaT 
i I i a 
PCO) CECH {| _0.7t015000 | 0.710 10000 | 
[Ampitude Range (3) SSS] SPR | OOSCSCdE:SC“‘OO:!OCO:OCSOOCOCSC“‘RSCSC‘“‘CSIN’“C(C*C*# 
[Resolution (oroadband) —sSs| gk | ~——o.o1—Ss«| ~~] SSCS SSC) SSSC«T=SC 
[Mechanical Shocktimits [gpk _{ 8000 | 5000 ~f 5000 § 2000  j{ 1000 
= zi 
| 181030 [181030 | 
P ConsiantCurentexciaion (6) | mA | 4 | 4 #4 4  § 4 +| 4 ~~ ~'| 
[Outputimpedance—S«Smm =O 1200~=CO«d:s=SsCS TOC’ 
[ OuiputBiasVonage S| SOC CTC Sto =| ~~~ tO TO =O Stoo ~OTS Sto CS 
[Discharge Time Consian— —=—=s«dt_second | 205 =| 205 Se aaa 
[Sensing Element ———=S=SSSC~*d=SCQuatz =| ~—eshear_ | ——irishear | 
P1082 coax 10-32 coax | 10-32 coax _| 
[weldinerm [_weldmmerm 
[Mountingthread «| Sse] AOTC CCS 
p Housing SSSSC*d CS mattoriel | manum [titanium of —sitanum sf sttanum [anu 
nore | aw wo oe | _-— a nn | -  | 
PTopconnecior | SCT SSSI ~—=Cstanddard —| Pssies2 [| ae1e42 
[ Groundisoiated(>10%0nms) | «dT ~SCSSTBTT ~—«|~—CSTBTA 
| Adhesve Mount | SS *dSCCASTBTT «| ~—CASSTBT 
SUPPLIEDACCESSORIES;(@°9) fC 
| “Mounting Sudo) | COUN COdTSCWNA Cd 
CEHONANACCESOORIEG:(}a0 | St Sit. | <a a” 
[Adhesive MountngBase =| —=~S=SS«d|~—~—COBOATS | ~-OBOATS 
[Magnec MountingBase «| df OBOASO =| ~—=COBOARO CS 
| Traxial Mounting Adaptor __—s«{_—~—=S—s«d|~__—OBOAT6 | —OBOATE 
CoN a (ee (ee ie nn Se oe | 
[Mating Cable Connectors «| ype =| SSG CUT SCCAC*S? 
[_RecommendedSiockCable (| series | 002 ~| ~~—~-oo2—=z 002 

Teo: Typical Sensitivity Deviation vs. Temperature 


. All models nave a sensitivity tolerance of +10%. 


Sensitivity determined with 4 mA constant current excitation. 

{Consult PCB for calibration at other current levels.) 

Amphtude range is for a 1.5 volt output full scale. 

In-house testing capabilities muted to -320 F (-196’C). Sensors 

nave deen successfully operated in quid helium enwronments. 

Within 5% of typical graph 

Operational range: 2:12 mA. Sensitivity drops 1% for each mA increase. 
See page 8 for a description ef the optional models. (Specitficabons may 
aitter slightly Consult factory before ordenng.) 

See cables/accessones section beginning on page 78 for a complete 
description of the supphed and optional accessories. 


Sensitivity Deviation 
«rel 100 Hz @ toom temp: 





Temperature ( F} 


9 NIST traceable calibration certiticate from 10 Hz to +5% frequency point 
supplied. Exact temperature coefficient at -320 °F {-196 °C) also provided. 
10. Supplied mounting stud is threaded 10-32 to 10-32. 


wl 


APPENDIX A4. ECU/ACTUATOR INTERCONNECT DIAGRAM 
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APPENDIX A5. GEOPHONE SPECIFICATIONS 
(Made by Halliburton Geophysical Services, Inc.) 


SPECIFICATIONS FOR SM-6 GEOPHONES 


Frequency 

Natural Frequency 

Frequency Tolerance 

Maximum tit angle for specified 
Fn (vertical units only) 

Typical spurious frequency 


Distortion 

Distortion with 0.7 in/s pp coil 
to case movement 

Distortion measured at 
frequency of 

Maximum tilt angle for distortion 
specification (vertical units only) 


Damping 

Open circuit damping for 
3750hmcoil . 

Shunt resistance used for 
damping calibration with 
375 ohm coil 

Damping with shunt for 
375 ohm coil 

Damping tolerance 


Sensitivity for Model ‘A’ Virals 
: Viin/s 

Sensitivity for Model ‘B’ B+ Vimis 
Viin/s 


RtBcFn for Model ‘A’ 
RtBcFn for Model'B' 


ee ee 
ee ee 
a 



















2 model guaranteed for one year Suspended mass Moos! ‘A’ 16.19 Key 
é i } Model “B' 11,19 " 2G 
a/ use excludes damage caused by high voltage and physical damage to case Teeltbaniee ee mesh ny ee perewicy ap 
ae . Operabng temperature range -40°C 1b + 100°C «= Open Crcut Cam 5 PLAC 
Sensor reserves the right to alter specifications to offer the best possible product. |-40°F t +212°F) Cod Current Serie a ne 5, 
Al parameters specibed xi 20°C Cod Ressstancs Ae tre + Ris} 
Shunt Resistance :Rs 
Tow Dameang & = Bo+ Bc 
Tot Remstance sa Re + Rs 
RLBCF9 Model “A = 10.8 Re 
Model 'B ~ 16.Rc 
Transducton constant 


Model “A’ G x: 0.0367 Ve Vinvsec. 25% 1.445 / Pc Virvsec, 25% 
Mosel 8 G = 0.0377 Ac Viln/sec. £5% 1.484 / Pc Vim/aec. +5% 


1a, 


OUTPUT Vil/s 





The SM-6 Is a digital grade, miniature geophone based on the popular SM-4 design. 

Engineered for consistent performance over a long life. 

Advanced computerized testing, manufacturing techniques and quality control are used in the production process 
to provide both the uniform parameters and the rugged qualities necessary in modern geophones. 


Features of the SM-6 

@ Available coil movement of 4 mm (0.16 in) which 
makes it suitable in a vertical mode and high tilt for 
lower natural frequencies and for industrial purposes 
where the force of one G is exceeded. SM-6 FREQUENCY RESPONSE CURVES 

@ In the honzontal mode the large coil excursion of the 


10 Hz model makes it suitable for shear wave SESARC OETECTOR RESPONSE CURVE 
detection with planting faults of up to 20 degrees ee eee 


without resorting to cumbersome leveling devices. 

e Symmetrical rotating dual coil construction minimizing 
the force on the spring arms. 

e@ Use of precious metals ensures optimum electrical 
contact and long operating life. 

@ The unbreakable and lightweight PE-4, PE-5 and HPE 
case provide for easy access to the element in the 
field. 

e@ Sensor's UV stabilized rubber sleeved gasket 
minimizes cable fatigue and has excellent watersealing 
Capabilities. 


OUTPUT V//S 


i 


OUTPUT V/m/s. 
OUTPUT W/i/S 


oe 
— 
er 
ae 


OUTPUT Vimis 


QUTPUT Virmis 


OUTPUT Vil/S 


COMI er 
SOs Ne 


atl. 


OUTPUT Vim/s. 
OUTPUT Vi/S 





2 2 ee 7 iT} 20 bv OND DY 10 200 Os : 
FREQUENCY “2 FREQLENCY y 


ALL SALES AND DELIVERIES ARE SUBJECT TO SENSOR STANDARD TERMS AND CONDITIONS INCLUDING THE GUARANTEE PROVISION 
STATED THEREIN. 
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OUTPUT Vim/s. 


APPENDIX A6. SPS390 SETUP AND MENU SETTINGS 


File Menu 
Open Configuration > stdSchan (configuration settings are listed below) 
Save Extended Record ~> used to save .xrc data files 


Setup Menu 

Acquisition 
Coupling > AC 
FS rms > user selectable, typically 1-5volts 
mV/EU > 100 
EU Name > Units 
db Ref > 0 
V Ref > 0 
Source > External 
Repetitive > as desired 
Source Channel Delay > -/00 samples 
Sampling Method > Internal 
Frequency Span > 4kHz 
RPM Readout > deselected 

Analysis 
Channel...Active > select all & channels 
Channel...Reference > select channels I & 2 
Block Size > 4096 samples, 1600 lines 
Overlap > 0% 
Process Weighting > Hanning 
Memory > 8192 samples 
Domain > Cross Properties 
Method > Linear 
Stop Criterion > Single 

Preferences 
DDE Server Capability > selected 
Display Buttons > selected 
Controls on Top > selected 
Overload Audio Tone > 500Hz (as desired) 
Shorthand Notation > deselected 
Peak Interpolation > deselected 
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Setup Menu(cont) 
Display 
Function > Time Trace 
Source > Live 
X-Axis Units > Seconds (Linear) 
Y-Axis Units > V (Linear) 
Active Reference > channel 2 
Tics > deselected 
Grids > deselected 
.Shorthand Notation > deselected 


Mode Menu 
Mode > Cross Properties 
Waterfall > deselected 
Total Input Samples > 65536 


* It is necessary to open a compressed time trace window for one of the reference 
channels under the “Window” Menu for proper data formatting during save operations. 
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APPENDIX A7. JOHN DEERE TRAILGATOR SPECIFICATIONS 





Specifications. 





Engine and electrical 4x2 
NARS ook i. oc0'c.0 6 « so eee > Legegegs ARE SPI oN noo So | erie Ine re Sc): 
NCS 25 aan ePera oo 6520 3 ACO CER)AS 22,6) cc cerita orcs = MNS o oiains cos ws «a «elon AGOVOIE OOS os 25 vsisuu' da castes ee eee 
CVMIIGEES 2s a. 5/0 5: o'ie Siete ccs 0:5 = semen rere EIeNe 1c ASME so wo o.oo thw wide oe = wd Gis ree Meena ny Uitte = Sera 6 Gara eer te rent 
WalVINg @ cornet os 3s <2 eee ws oe eee DOVECNGaU MaIVES can ae fete ka's s Sein ea lew cdc ae ere OMEMIGAC WANVES ceecscaccen ors oo eee ee 
Displacement? 22 estes a ds oe eee Ae SCUIN (LOO ICC) anvines.c css -siettewe nuns eee eae BO led CUMING (ONE CC) o0.iea.20 3 eee 
Horsepower (SAE J1940)................ URS as ee mee rcmrrere 1a - cacch a eh Si OES ea 
Maximum tOrque.... 2... 00s. ccaccccreces MAE) Derk Cae, UDI 25 x os a'o'eslee seis se ses Ce ee 32.5 Ib-ff @ 2.200 1rDM. ... 2. ccc ccnswcscns ee 
DUG 65a ae mae ae erating irene Ga BaClnve CIGCN ANG ns <a <b. cate es cicwor a oem horeraeee Wieansistor Controlled. 222... 5 .casces vs ce ate 
PPADELGaUOM bs yee cise s.s sa eideleie eas cave cures RUlpres tle s crc e acon s28 8 dss dn weae os opie wis so tee IPO DROSS UNG samererney: oc. '.-o 2 ee 
LUTE TUT Se ene oe a cee Replaceable {optional}. ...........605.500s0es00ce.. | meplaoeable (Standard). ......:...<.4¢-)25se 
Rpm, idle (no load)............. eee UA Ae a rrr irs) | (||| 7c) Gamat Gry RR ee Go 
RDM fast CMONOAG) - csc. x we cee es ene. SUM meee ieee ea ies aia ae Coe vin ww «te ee PO Uer.s s On cares oe cae ee eee ee 
GOGINOISYSIGIN <cn.< vies seen «enn ETS Et rar eee PRTG et eS EN 2 A Nias, 3. «nse rn ance lw We'd Savoie eee 
1s \ifitcel [9 fal a Ory, replaceable, single element, with remote intake ....... Dry, replaceable, single element, with remote intake 
LET a a PS rarer ee SPalicaire sting: eect cs cs. caus. eee SPAM AaMeStNGerscs soe eres oe «os xls P 
BPA Peete ie ores atte cia teiaiely si¢.s'svo aya che eee 32-amp/hr, 295 cold cranking amps ..............000. 32-amp/hr, 295 cold cranking amps ........... 
PULCRALOlitctt <2) c:0:5 0 2 os ere cee kivis a scine ee TS-aMP TEQUiated 6c. occ so. co cee cde wce oe adiess ee eam MC OUIALEO Mer err. 5 oc. caveats eee 
PP AUGMISM er eck cece cok kee c cas TWO137 O-WallHdlOGeN «2.2... oig.anunswcsda ssa ueenee Two 37.5-watt Halogen ..................0.. 
Fuel system capacity................... Did We OnUaCe lier Ot «nade Wan ee esate meen io Us Gale (ZO. + c.0ais os teneinles seer 
Transmission 
FUT meme epee itty cie-c)» 1.0 cee sas « eueuaitpsieece’ Continuously variable transmission (CVT) ............... Continuously variable transmission (CVT)....... 
Bitterentiailock ta... 2. .erotnes e « Standard; hand operated .......... PT Ne Standard; hand operated .................... 
DPLETIAI eer eas aiwinx os seo ees BE: 2 Nae Nt, Sos ae Mh eaten .... #No. 50 industrial roller chain ................. 
GWUNS ee Ga ncaa cdascies we oe ae O=TS-OIMPNO-2a AMUN)... os ened ekd ss es cose e wel se Sle ee (250 1)/)) i re 
HiPABL SANE ER Seti 5s sy Pade aee os cea Fully enclosedvonetiat. os .c46 saci enc scree ae eeeden Fully enclosed; oil bath...................... 
Bede SblechOnes:. 22. %..0> os shee st ws FOlWald: NeUtral MEVENSe .ccuesd os wale ete aula ue FOrWard MeUtalameVverse .. 22.0632 sane ew ciaw's 
Overall reduction ratio: ; 

WOW ee cet cs su creineSe cies ae Sees Neer ere een sa oon Nh eis ct ee ee TES <A oe Sn en enn ROMER Inn ey Ni. 

Ble Serres toes yr pee ce arene ROIS: lene eee cco. kc nue taees eae cae ee A Gio te ss cs os s rae eee ee omer 
BMS oceans seein mares Wet disk in transaxle ......... cece eee eee cece cece es fp Wet Gisk in transaxle 2.0... eee eee 
SS) CLL LED Epson oan En en Rena ser Rack and pinion with Ackerman-type geometry.......... Rack and pinion with Ackerman-type geometry .. . 
Ground pressure {max.) 

With 200-lb. operator only ............. SPOIDSH (SRO GINS) cc anes aslvin Gwie oar eee see BONS CONA) 6 sss nse oc. Seed desu 

Fully loaded vehicle................... TiS PSI SO/ KOU CIN) «2:4. s < ga cinnwcies soatev eles ce sane HON DSIEQUIRO CMS) cis. o as cs so sco du cine ees 
Dimensions 
Length (without bumper) ................ MOO: Suet Zoo CM). cre cs te ee oe wer ov deg sar OAS INIZ6610 CM ns hs ese ode deena oe 
Length (with bumper)................... IOLG In e56em). 02s. oe ees oe cee O in. SC) ee co eee 
CUT RRs 3 one COMIN (ozcIGIN) <2. ewe neers ccs oe au cee eee .0 in. BSTC IM c-concee have eo oeie hee cate 
FrOMi-nede COMICIS 5. sn es ccc ccc seeceees DOO IMURECeGM) cc 4 ocic vcdicin ac css crceles oss Sos faces OITA oc swans oc Sle os oe otete inet 
Rear-WeaG CENTS . .se.s- ence cee se scene MINS CIGMD ie. ..c ccecsleiale wicicis treed diwce aves eis sae ae MOIR IZe GIN) 2 <5 6 ssn = gic ee coe eee eee 
HEIGNULOVETAIN) 5... cece coe cs aetewenect MOPOIIM MMOL? (GIN) ct uiciare ttelo s Berns ceeds esa ual 6 in. BITCINS Sree s «tee Rance eta Rien. ale 
PCC DAS Oe ati otic orate eae walside eee GS: Sime (UF6.6.CM).. . scare ck cc peace bas oc os ces cece eHeeO INmteoe oe OM) MONUIeal. «cscs. bce ees cae 
Weight (including fuel/fluids) ............. B/E GSS OM) cco es ie eee as. bs bac abcde aides ois AMMRUOU DANA DOS IKO) So osc vw cane eeinnea nauee 
Ground clearance: 

Under transaxle............ ee ier Grama ROMM ire nc-c eases ncaa soy ese ee eee Cee ALON) siscectvs sco tcaroncare a ange ene 
Near {OOt PAVONN 4... occ ees ss os eee es S Sin CacM eo eee ey ects ut ease eee ae Bre OWEN) lire «69 Geass Gols e werstee eee we 

SPalinG Capacliyiee «acess a. s 6 oes eos 6 oo Pa srg ich a a ee. CO cn a ee ee a 
STATON 6 oc 2 9 AS ee eee Standard medium Dack. . <6... 0c. ce cece ces scecce ssc cmmeratessional nigh back... ....2 6. ecccscocs cove 
urn cleatance Cc’. 0... coerce n. +s oeese Ze VRC MU Me elaceicrs «vo ave ecuelts Coie ses eels eee oat 30) 212A SN) le i ae 
MOWINGICADACIIY ae-4s 5's. sa ass vs caesar MUO MAUS KO) Gory neces ecercten 6. see's sence : EL) a ar ee eee eae 
RavloadiedQacny wasn... sseensaasehees SOO IE AOS KO) iss cai: sieilas er ars s oie. 0eks «0 <a re em COUN OMe KO): 5 6.4 5-5. ce Coo te oe 
Cargo Box 
Maternal ose eck war ees ne eee GI MICIS TEEN oo oe Croc totes ele sc Sais iea a oc so sc ree REE OEN We iv n So 5c os s5e dew eves ee ame mts 
Capacity: 

VOIUIM Ecce ca eer tae sven nee AOMBCHMICZ SING) 5. 0sa Seo cine s bees seve se cecenme EC OURMOUAZ NY) 5 2-2 + ge we es neoe nels 

WIQNT rey cs oes ve ee SUMMER ere co ace hacseeree eee, ale a vies 9 wa ee MCRD) cs nn eos we eas eee eee 


"Includes 200-pound operator, 200-pound passenger, and maximum box capacity (Specifications and design subject to change without notice.) 


me] 








APPENDIX B. UNDERWATER ACTUATOR CABLES AND CONNECTORS 


This appendix contains the specifications and wiring diagrams for the underwater 


actuator cables and connectors that were used during this research. 


Appendix B1 — Belden 9775 cable specifications 

Appendix B2 — Seacon rubber-molded connector specifications 
Appendix B3 — Seacon bulkhead connector specifications 
Appendix B4 — Underwater cable wiring diagram 
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APPENDIX B1. ACTUATOR UNDERWATER CABLE SPECIFICATIONS 


— ee 
Paired [Belden 
Individually Shielded 





Audio and Computer Cables 


Standard Nomnat No-runal 
Lengins Nonunat 0 C.R. 0.0. MOMs, Cavacitance 


Trace No. No : 
iuL NEC Tyne of 
CSACert. Paws Lbs 


Descmption 
ur AWM Style ° & es es 






ft m 2. Cord Sthieta inch mm ails. Prop : = ed ae 
18 Gage Product Description 
Stranded Conductors (19x30) Tinned copper, polyethyiéne insulated, twisted pairs. Each pair individually shielded with 
Polyethylene Insulated Beldtoil® alumninum-potyester shield and 20 AWG stranded tinned copper drain wire. Overa!| 


chrome PVC jacket. Color code chart No. 3, Technical Information Section. 















Beldfoit 
100% Shield 
Coverage 
WA 2919 
30V 80°C 


See Attenuation, 
Rise Time and 
Bit Rate data for 
this series on 
page 99. 


Belden 9775 
Cable used 


* Capactance between conductors 
102 “° Capacitance between one conductor and other conductors connected to sheld. 
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APPENDIX B2. UNDERWATER RUBBER-MOLDED CONNECTORS 


RMK-FS AND ASSOCIATED PARTS 
STANDARD RUBBER MOLDED 


IN-LINE SERIES 


ey ee 
Connectors 

Material: Special Neoprene, Molded 

Contacts: Copper Alloy, Hard Gold Plated per MIL-G-45204 

Voltage Rating: 600 VDC-contact factory for higher voltage requirements 















Pressure Rating, Mated: up to 20,000 psi ? - 
moex 
| A. = = 
| cick we 
He 
RMK-FS — 


Co eo ret ek > a ot 


ae 
Ls 


4 
ea 


1 


RMK-FS R/A 


the, 


°o ee 
ee ce) 


Rae pal RMAK-1-FSAMP Lite d RMK-2-FSMP  RMKSFSP | ORMK-PSAIP  ORAIKS-FSAAP ORMEFSP | RMK-7.FSMP AWS 
14823 1810 1410 41 #10 1410 vara 14 


egw 
(Care: —_— (Cabis: 87150) ye 2014 gai aes 6016 7416 
= 2280) (Cable: 18/250) (Cable 18/3S0) (Cable: 184S0) [Cabe: 10550) (A FESO; canta 187750) (Cable: 18/850 


LOCKING SLEEVES DUMMY¥rCONNECTORS 
RIMAK-FS to a bulkhead connector. Must specify cable OD. 

K-MLS-D (AVAILABLE FOR RMK-MP) RMK-FSD MATES TO RMK-MP 

K-FLS-D (AVAILABLE FOR RMK-FS, RMK-FS-R/A, WOC & BCL’S) RMK-MPD MATES TO RMK-FS & RMK-FS-R/A 


Matenal: Special Neoprene, Molded 

Pressure Rating, Mated: 

RMK-FSD (FOR HIGH OR LOW PRESSURE, NO CONTACT REQUIRED) 
RMK-MPD-LP (LOW PRESSURE < 4,000 PS!, NO CONTACT REQUIRED) 
RMK-MPD-HP (HIGH PRESSURE > 4,000 PS! CONTACT IS REQUIRED) 


WOC CONNECTORS: 


Note: WOC stands for “without cable.” The connectors are supplied as shown below for installation in the field 


and/or to non-neoprene jacketed cables. Bonding can be achieved by molding to the connectors shown using 
polyurethane with appropriate primers. Inexpensive plastic syringes can be used as molds. 





Part Names: RMK-FS WOC, RMK-MP WOC 

Material: Special Neoprene, Molded 

Contacts: Brass, Hard Gold Plated per MIL-G-45204 

Voltage Rating: 600 VDC-~-contact factory for higher voltage requirements 
Pressure Rating, Mated: up to 20,000 psi 
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APPENDIX B3. UNDERWATER BULKHEAD CONNECTORS 


RMK-FS AND ASSOCIATED PARTS 
MATEABLE BULKHEAD CONNECTORS 


BULKHEAD CONNECTORS 
(MATES TO RMK-FS & RMK-FS-R/A)} 


Part Name: RMK-BCL 

Maternal: 316 Stainless St.; Intemal Insulation: Special Neoprene 
Contacts: Copper Alloy, Hard Gold Plated per MIL-G-45204 
Pressure rating, Mated: up to 20,000 psi 


Pan Name: XSK-BCL 

Material: Glass Reinforced Epoxy 

Contacts: Copper Alloy, Hard Gold Plated per MIL-G-45204 
Pressure Rating, Mated: up to 20,000 psi 


VISUAL INDEX 
(N-UNE WITH CONTACT #1) 








RAMK-6-BCL 










RMK-2-BCL RMK-S-BCL 
XSK-2-BCL XSK-5-BCL XSK-6-BCL 

1414 1010 1014 

1410 4014 § 016 
Part Name: XSK-DEBCL 
Material: Glass Reinforced Epoxy 
Contacts: Copper Alloy, Hard Gold Plated per MIL-G-45204 
Pressure Rating, Mated: up to 20,000 psi frorn front side. 
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XSK-DEBCL 


Description: Special insulating neoprene is 
permanently molded to the meta! body. These 
types of connectors are used where physical 
abuse would be too great for our glass reinforced 
(GRE) epoxy series. 


Description: A non-metallic connector 
homogeneously molded of glass reinforced epoxy 
{GRE). One of the oldest and most reliable 
connectors in our product line. For the same 
connector body, but with female sockets refer to 
the VSK (page 50-31). 


Mounting Torque VISUAL 
130 inch Ibs. INDEX 
é ¢ g 
RMK-1/0-BCL RIK: 1-BCL PAIK-6/8-3CL 
XSK-1/0-8CL XSK-1-BCL XSK4/8-BCL 
1410 183 106 


198 





RMK-7-BCL . 
wees ae hint 
4 re 7 316 11 #16 
Vv 
RMK-9-BCL 
Connector 
Used 
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APPENDIX B4. WIRING DIAGRAM FOR UNDERWATER ACTUATOR 


CABLES AND CONNECTORS 
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APPENDIX C. EXPERIMENTAL SETUPS 


This appendix contains the experimental setups for selected tests presented in this 


thesis. 


Appendix Cl — Ambient Noise Analysis 
Appendix C2 — Vertically Oriented Seismic Source Setup 


Appendix C3 — Rayleigh Wave Seismic Source Setup 
Appendix C4 — Frequency Optimization and Wavespeed Measurement Setup (Vertical) 


Appendix C5 — Test setup for helium gas tank target tests 
Appendix C6 — Test setup for helium gas tank and gunpowder keg target measurements 


(increasing mass) 
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APPENDIX Cl. AMBIENT NOISE ANALYSIS 


This appendix contairis the experimental setup for the ambient acoustical noise 


analysis conducted on September 4, 1998. 


Comments: Geophone gain = 40db 
Sunny, low tide, 3-4ft waves. 
Unfiltered signals were recorded as shown below 


> e 
a5 Waterline 
RRR VV 


8 b, 


a 110 





86 


APPENDIX C2. TEST SETUP FOR VERTICALLY ORIENTED SEISMIC SOURCES 


This appendix contains the experimental setup for vertically oriented seismo-acoustic 


sources on October 20, 1998. 


Source #1 Source #2 

Orientation: Vertical Orientation: Vertical 
Voltage: 10V (20Vpk-pk) Voltage: 10V (20Vpk-pk) 
Frequency: 90Hz Frequency: 90Hz 

Cycles: 1 Cycles: I 

Geophone #1 Geophone #2 

Filter:. High Pass 40Hz Filter: High Pass 40Hz 

Gain: 40db Gain: 40db 


Comments: Sunny, low tide, 3-4ft waves. 


60 | Waterline 
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APPENDIX C3. RAYLEIGH WAVE SEISMIC SOURCE SETUP (PHASING) 


This appendix contains the experimental setup for Rayleigh wave excitation on 


October 2, 1998. 


Source #1 Source #2 
Orientation: Horizontal Orientation: Vertical 
Voltage: 6V (12Vpk-pk) Voltage: 6V (12Vpk-pk) 
Frequency: 100Hz Frequency: 100Hz 
Cycles: iD Cycles: 2 
Phase: O- Phase: -90° 

Geophone #1 Geophone #2 
Filter: High Pass 40Hz Filter: High Pass 40Hz 
Gain: 40db Gain: 40db 


Comments: Sunny, low-medium tide, 2-3ft waves. 









45 Waterline 
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APPENDIX C4. FREQUENCY OPTIMIZATION AND WAVESPEED SETUP 


This appendix contains the experimental setup for determining the optimal drive 


frequency and its associated wavespeed on October 20, 1998. 


Source #1 Source #2 
Orientation: Vertical Orientation: Vertical 
Voltage: LOV (20Vpk-pk) Voltage: 1OV (20Vpk-pk) 
Frequency: 90Hz © Frequency: 90Hz 
Cycles: St Cycles: ] 

Geophone #1 Geophone #2 
Filter: High Pass 40Hz Filter: High Pass 40Hz 
Gain: 40db | Gain: 40db 


Comments: Sunny, low tide, 3-4ft waves. 


| 60 Waterline 
ARR RR TTS 


G2 


( }+—> 





Not to Scale 


89 


APPENDIX C5. TEST SETUP FOR HELIUM GAS TANK TARGET TESTS 


This appendix contains the experimental setup for target detection with increasing 


mass tests conducted on October 23, 1998. 


Source #1 Source #2 
Orientation: Vertical Orientation: Vertical 
Voltage: 10V (20Vpk-pk) Voltage: 10V (20Vpk-pk) 
Frequency, 80H Frequency: 80Hz 
Cycles: i Cycles: ] 

Geophone #1 Geophone #2 
Filter: High Pass 40Hz | Filter: High Pass 40Hz 
Gain: 40db Gain: 40db 


Comments: Overcast, med-high tide, 4-Sft waves. 
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APPENDIX C6. TEST SETUP FOR HELIUM GAS TANK AND GUNPOWDER KEG 
TARGET TESTS (INCREASING MASS) 


This appendix contains the experimental setup for target detection with increasing 


mass tests conducted on November 6, 1998 and November 10, 1998. 


Source #1 Source #2 
Orientation: Vertical Orientation: Vertical 
Voltage: 1OV (20Vpk-pk) Voltage: 1OV (20Vpk-pk) 
Frequency: 80Hz Frequency: 80Hz 
Cycles: J “Gyles: 

Geophone #1 Geophone #2 
Filter: High Pass 40Hz Filter: High Pass 40Hz 
Gain: 40db Gain: 40db 


Comments: Overcast, med-high tide, 4-5ft waves. 


Helium Tank (Nov 6") 


Gunpowder Keg (Nov 10") 
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APPENDIX D. COMPUTER PROGRAMS AND SUBROUTINES 


This appendix contains the Matlab subroutines and programs that were used 
during this research for signal processing and data analysis of the recorded beach data. 
To obtain a copy of the subroutines contact either Professor Thomas G. Muir, Code 
PH/Mt, or Professor Steven R. Baker, Code PH/Ba at the Naval Postgraduate School, 
Monterey, California 93943-5000 


Appendix D1 — Attenuate.m 
Appendix D2 — Avgsum. . 
Appendix D3 — ChanI.m 
Appendix D4 — Chan8.m 
Appendix D5 — Coherence.m 
Appendix D6 — Convert.m 
Appendix D7 —Correlate.m . 
Appendix D8 — Displace.m 
Appendix D9 — Filtr.m 
Appendix D10 — Frequency.m 
Appendix D11 — Hankel.m 
Appendix D12 — Hilbrt.m 
Appendix D13 — Hodogram.m 
Appendix D14 — Modification.m 
Appendix D15 — Patch8.m 
Appendix D16 — Powersd.m 
Appendix D17 — Rms8.m 
Appendix D18 — Subhilbrt.m 
Appendix D19 — Subtract.m 
Appendix D20 — Velocity.m 
Appendix D21 — Wavespd.m 
Appendix D22 — Xzplot.m 
Appendix D23 — Theory_Displace.m 
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